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We present a complete classification of the electron-electron interaction in chaotic quantum dots 
based on expansion in inverse powers of 1/M, the number of the electron states in the Thouless 
window, M ~ IzfR- This classification is quite universal and extends and enlarges the universal non 
interacting RMT statistical ensembles. We show that existing Coulomb blockade peak spacing data 
for B = and B 7^ is described quite accurately by the interacting GSE and by its extension to 
B 7^ 0. The bimodal structure existing in the interacting GUE case is completely washed out by 
the combined effect of the spin orbit, pairing and higher order residual interactions. 



Measurement of the fluctuations of the Coulomb block- 
ade (CB) peak spacings as well as the mesoscopic spin 
effects provide an excellent probe of the properties of 
interaction effects in disordered quantum dots (QD). Ba- 
sic experiments where performed |IL^| in ballistic, well 
isolated QD of irregular shape formed in 2d GaAs het- 
erostructures. The QD sizes were R 3> 1/kp, where kp 
is the electron Fermi wave vector. 

The main general theme underlying this problem is an 
interplay between chaos and interactions in quantum me- 
chanical motion of electrons in a restricted geometry. In 
a recent work || , cf. , also Q , on mesoscopic spin effects 
a proposal appeared for a universal Hamiltonian which 
controls the main physics of interactions in a chaotic QD 
in the extreme limit kpR S> 1. In || statistical fluctu- 
ations around this limit were included in the framework 
of the Hartree-Fock (HF) method. Our main goal here 
is to show that these fluctuations and indeed the entire 
interacting Hamiltonian can be represented and classified 
in a unifying universal scheme, cf. Eqs. below, 
which follows and extends the universal symmetry classes 
of the Wigner-Dyson statistical theory for non interact- 
ing electrons. Our second goal is to apply this theory to 
the problem of fluctuations of the CB peak spacings. Al- 
though this has received much attention recently , 
a consistent description is still lacking. We will show that 
the universal interacting GSE Hamiltonian and its exten- 
sion (we call it GUSE) to the non zero magnetic field in 
our scheme accounts quite well for the experimental dis- 
tributions, cf. Fig. 1 below. The GSE choice matches 
perfectly the recently discovered strong spin-orbit (SO) 
effects in GaAs QD, [| 

The Hamiltonian of an interacting QD consists of one 
and two body parts, H = H n + H int . 



Ha— ^ Hao,b<j'a\cr a bcj' — ^ EaaJ^cw, (1) 
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Here indices a, b denote space orbitals while a, a' are the 
spin indices. With a possible SO interaction H a<T ^ a i is in 
general not diagonal in the spin indices. We use a, v to 
numerate the eigen states of H aa ^ a i ; e a are its eigen ener- 
gies with a the orbital and v = ±4 - the spin or Kramers 
index in the presence of SO. In irregular QD the one elec- 
tron Hamiltonian in the Thouless window of states can 
be described by a random matrix theory (RMT), [[To] fl2| . 
We will denote by M the rank of H ar7 i„j> . The statistics of 
H a <j,bcr> depends on the symmetry of the problem classi- 
fied by standard RMT ensembles - GOE, GUE and GSE. 
The correlators of the eigen functions of R a a,ba' depend 
on the ensemble. For GOE they are 

< ip* au (r,a)il)^(r',a') >=< </w(r, <x)^'(r', a') > (3) 
=< ^>*„(r, a)ipp v ,(v' >= 5 ya 8^ al 5 aP K{v,v'). 

We find it convenient to work in the coordinate-spin rep- 
resentation (r, a). The function K(r,r') is 



M 

K ( r > r ') = ^EE ^( r > °)^(r', a) (4) 
~A- 1 J (fc F |r-r'|). 



Here A is the area of QD and Jq(x) is the zero order 
Bessel function giving an approximate quasiclassical ex- 
pression for K (r, r'), For the GUE the correlators of 
the type < iptp >, and < ip*ip* > are zero while < ip*ip > 
is the same as in GOE. For the GSE symmetry one has 



< ip* au {T,o-)^p vl (v' ,a') >= -5 1/v >dv0>d a pK(T,T'), (5) 

< '4> a u{r,o)ipfiv{r\a') >= -AW'AW<W-K(r,r'), 

k is 2 x 2 time inversion matrix for spin ^ systems p0[ ] . 

It is convenient to consider the interaction part of the 
Hamiltonian Hi nt in the basis of the eigen functions i/w 
of the one electron Hamiltonian 
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( r > ct)-0^ 2 (r', o')U (r, r') W.3 0, (r', a')- 

Here (r, r') is the the screened electron - electron inter- 
action in QD. If SO interaction is absent the space and 
spin coordinates r and tr are separated. 

It is possible to represent the interaction (^|) as a sum 
of parts of different order in a small parameter 1/M. An 
essential step in this direction was made in M. Here we 
shall present the complete classification and investigate 
some of its consequences. We will only present our main 
results deferring the detailed derivations to Ref. [p|. It 
will be sufficient here to assume that higher correlators 
of V>'s obey the rules of the Gaussian statistics. The role 
of the non Gaussian corrections will be discussed in |L5] . 
We use a cluster decomposition of the matrix elements (H) 
as fourth order polynomial functions of the random wave 
functions ip. Inserting this into Hi nt we find that it con- 



sists of three groups of terms H int 
For the GOE these terms are 
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Hl°l = -V C N(N - 1) - JS 2 + Pf tT, 

= Nj2<p(at • o/j) + S • <p(al<ra p ) 

a/3 a/3 

P 



(7) 



ap a/3 
(2) 1 ^ — ~ t t 

H int = y 2^ Vap-jsal^a^as^a^. 
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Here we denoted by dot the scalar product in the spin 
variables, TV = J2 a ( a a ' a a) is the operator of the to- 
tal number of electrons, S = i J2 a ( a a cra a) 1S the op- 
erator of the total spin, er are the Pauli matrices and 
T = J2 a a aia a ^ is the total pairing operator. The de- 
composition (Q) is an identity which is useful since as we 
will show below and in [|l5| it allows to classify the groups 
of terms , H^l an d H^l by their degree of smallness 
with respect to 1/M. Namely, apart of the capacitance 
term V c ~ MA (A is mean level spacing) and up to 
logarithmic corrections the matrix elements in H^ t are 
~ A/(Af) n / 2 . The properties of u and V are discussed 
below,cf. also Ref. |.|. 

For GUE the expression above remains valid except 
that the terms containing the pairing operators T and 
are absent. For GSE the expression (pi) becomes 



H%l = -V C N(N-1) + Pflf, 

OLUpv' olvQv' 
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H {2) = - 
int 2 
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Now the terms containing the spin operators disappear 
and the spin index is replaced by the Kramers degener- 
acy index which we denote by v. The matrices u n,p and 
V depend on v. 

Introduction of a perpendicular magnetic field removes 
the Kramers degeneracy and changes the statistics of the 
single particle Hamiltonian into GUE(2M). But this is 
not the interacting GUE obtained from (0). The second 
correlator in (||) vanishes and consequently the terms con- 
taining the pairing operators T, in (Q) disappear. Also 
there is no need anymore for the Kramers spinor indices 
in u n and V. Thus one obtains a different 1/M expansion 
which we term GUSE (unitary arising from simplectic) 



H 



(GUSE) 
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-V C N(N-1) + N u n a0 aiap+ (9) 



a, /3=1 



j 2M 

;7<5 a a a /5 



a, /3, 7, 5=1 



Such an ensemble was first discussed in || . 

The lowest order H^ t in (fjj) and its universality was 
fully discussed in g. The \V C N{N - 1) term is at the 

basis of the simplest Coulomb blockade theory while JS 2 
appeared in relation to mesoscopic spin fluctuations, [H. 
Here we will discuss higher order terms and will then ex- 
plore their effect and interplay with the T'T term. We 
note, that u and V are such that < u^g > = 
< Va/3~f8 >— for i = n, s,p, and the average of the sec- 
ond functional derivatives of V with respect to i\> or ip* 
is equal to zero. Up to corrections of the order 1/M the 
matrices u l aB , u l avBv , for i = n, s,p, are Gaussian random 
variables with zero average. For GOE one finds, P,|l5| , 



l ap a a'P 
v2 



(10) 
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Here Cu> — Cin are 6 dimensionless constants which de- 
pend on the average geometry of as well as on the details 
of the electron screening in QD (cf. , below) . On the ba- 
sis of the correlators (^) one can find similar averages for 
GUE, GSE and GUSE. For GUE the correlators of the 
residual interaction V in (H) are, (5j|ll| 

Tr* f> A 2 In M 
< V aP: 1 8 V a'(3'-r,'8' >= X 

xD{5 aa * 5pp> <5 7 y 5ss> + 5 a p>5p a i8 lS ,5s 1 f + (11) 
+8 aa >Spp>8~ / 8'88y + Sap'Spa'S^'Sss'} + 0(M~ ) 

where D is again a dimensionless largely universal con- 
stant (cf., below). One can easily write corresponding 
expressions for other ensembles. 
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In order to proceed it is convenient to adopt the fol- 
lowing decomposition of the basic screened e-e interaction 
U(v,v) in (§), cf., § 

U(r, r') = V c + u sur (r) + u sur (r') + V(r, r'). (12) 

Here V c is the constant capacitance part, cf., 
u sur (r) is the surface part of the potential caused by the 
screening charges which are on the surface and V(r, r') 
is the screened bulk e-e interaction. One can express the 
constants entering in the expressions in terms of 

u sur (r) and V(r, r'). For GOE one finds (expressions for 
other ensembles are very similar), Jig] 



J = P = h 



A 71 - 2 J K n (r, r')V(r, r')drdr'; 



(13) 



l a p = / w 2 (r,r')V'*(r)V'/3(r')drdr',i = n,s,p; 



Here u n (r,r') = 5(r - r')(u sur (r) - u) - (l/2)u s (r, r'), 
u p (r, r') = it s (r, r') , u = A" 1 J u sur {v)dv and u s (r, r') = 
K(r, r')V(r, r') — JS(r — r'). The matrix V a p- yl s is given 
by (^]) if we substitute U(r,r') U(r,r') and extract 
the irreducible part. We also have 



D = 



M 
A 

M 2 



c «' = -X2 I K ^ 3 )^( 2 > 4 K(u 2 V ( 3 < 4 ) dfi < ( 14 ) 



A 2 InM 



K 2 (l, 3)K 2 (2, 4)V(1, 2)V(3, 4)dfi, 



where F(l,2) = V(ri,r2), etc., and dfi = dri<ir2<ir3<ir4. 

We now turn to the problem of the fluctuations of the 
spacings between Coulomb blockade peaks. We focus on 
the experiments in 2D GaAs dots, Ref. Q. An important 
observation made in Q was that the basic non interact- 
ing Hamiltonian for such dots must include a strong SO 
interaction, the so called Rashba term |l4j, «soP' [s x n], 
where p and s are the momentum and spin operators, n 
is the vector of the normal to the QD plane. The strength 
of this term in a typical GaAs/GaAlAs heterostructure 
is aso — 2.5 • 10~ 7 mev cm/Ti [jL4| . The corresponding 
energy scale ~ 0.3 mev is ^> A. Thus it is appropriate 
to use the GSE ensemble for the random single electron 
Hamiltonian and the expression (||) for the interaction. 
The pairing term T^T unlike other zeroth order terms 
does not commute with the random single electron part 
and should therefore increase the effect of fluctuations. 
Experiments in |j| included also the situation with an 
applied weak perpendicular magnetic field. This corre- 
sponds to the GUSE Eq. (§). 

For the QD parameters we use A = 2h 2 /(m*R 2 ), 
Ep = irh 2 /(m*r 2 c ), r sc ~ n c 1 ^ 2 , where m* is the effec- 
tive mass, n c is the electron concentration in a QD. The 
Thouless energy is Eph = h~/ T bai — V2~7rh 2 / '(m* r SC R) . 
From Eq. (S) it follows that the rank of RMT is 
M ~ TTRk F /2 so that M ~ tt z / 2 R/ (2r sc ). The constants 
in the interacting part of the GSE Hamiltonian (p|) are 



completely determined by u sur (r), V(r, r') and K(r,r'), 
Eq.(|). We take u sur (r) = -(e 2 /e*KR 2 ){l - r 2 /R 2 )~ 1 / 2 
which is appropriate for a 2D disc of radius R in the 
limit kR 3> 1 where k — . For a disk shape one gets 
an estimate C = 2Re*/ir, V c = e 2 /C = MA/(2 v / 2^). 
The screened interaction V(r,r') must behave as e 2 /er 
for Kr -C 1 and as e 2 /en 2 r 3 for kt 3> 1, r = |r — r'|. The 
constants P, C and D in Eqs.(0|jl^) can be expressed pl| 
in terms of the integrals /„ (|13|), for n = 0, 1, 2. They are 
sensitive to the intermediate range behavior of V. We 
estimated them as I = 1.5A, h = 0.51A, I 2 = 0.37A. 

We treated the last term in ((t||||]) in the Hartree-Fock 
approximation. This and the term NY2 u ^v3v' a av a /3i'' 
in (|h lead to a modified single particle part of H 



ttHF 
H ab 



Nu 



ab 



E 



cd 



VabcdPdc 



where p is one particle density matrix, V^ cd 



Vacdb — 

V ac bd and we omitted the spin indices. One can show, 
fUsf , that statistical properties of the HF eigen values 
and of the corresponding eigen functions are practically 
the same as in the original RMT. The only noticeable 
effects appear when the particle-hole energy differences 
are of order A/M, |@. 

We have calculated A 2 (N) = E N+1 + En-i — 2E N 
in the GSE and GUSE cases and compared with the ex- 
perimental data of S.R. Patel et al., [Q. Here En is the 
ground state energy of QD with TV electrons. The results 
are shown in Fig.l. Our calculations in obtaining these 
distributions were kept at a very simple level. The GUSE 
was the simplest case since it did not have non trivial in- 
teraction terms in the leading 1/M order, Eq. (||). We 
used the HF expressions for E(N) and obtained 



\GUSE 
^2 



(N) = V C + e N+1 -e N + Wjv+i,iv+i + u n,n (15) 




.5 -1 -0.5 0.5 1 1.5 

FIG. 1. v = (A2— < A2 >)/A - normalized peak spacings 
for B = and S / 0. Histograms are experimental data, 
solid lines - predictions of the interacting GSE (GUSE) for 
B = (B 0) 



3 



We used RMT GUE distribution for ejv+i — ejv and the 



Gaussian distribution for u 



N+1,N+V U N,N 



with the co- 



variance C nn = 0.069. This value as well as C np and C pp 
below were obtained using a reasonable parametrization 
of the screened e-e interaction, [ ^5[ . 

For GSE the calculations required a proper treatment 
of the pairing PT'T interaction appearing in the leading 
order in (||). This problem has an exact solution Jl6| . 
However we used a simple approximation which we felt 
was satisfactory. For even N we minimized this term 
in the subspace of two HF solutions with adjacent filled 
and empty Kramers pairs. We then formed an expec- 
tation of H with the resulting wave function. For the 
odd N the effect of the pairing is much simpler and the 
expectation with lowest energy HF wave function was 
sufficient. The resulting expressions are too cumbersome 
to record, cf., @. We used them with the RMT GSE 
statistics, P = 0.37A, the covariance C nn as in GUSE 
and C pp = 0.023, C np = 0. 
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FIG. 2. Normalized peak spacings for the interacting GUE 

As one can see in Fig. 1 there no sign of the bimodal 
structure in the GSE and GUSE distributions. The rea- 
son for this in GUSE is perfectly obvious, cf., Eq. Jl5| ) 
- the spin degeneracy which is responsible for the bi- 
modal structure in simple models of the CB is completely 
washed out by the combined effect of the SO interaction 
and the magnetic field. In the interacting GSE the non 
commutativity of the paring term PT^T with the single 
particle part causes rather strong fluctuations relative to 
the RMT already in the lowest order. To appreciate this 
effect it is instructive to compare the leading GSE inter- 
action terms, upper line in (H) with those of GUE, the 
upper line with P = in (j?]) . The commuting spin inter- 
action JS 2 does not change the basic RMT fluctuations 



but simply cuts and shifts different spin parts creating 
sharp structures. As seen in Fig. 2 these structures are 
washed out only partially by higher order terms. We con- 
clude by observing that our results, Fig.l, fit quite poorly 
the tails of the spacing distributions. It is not clear to us 
if this is a consequence of our approximations in calcu- 
lating Aa(iV) or because of more fundamental reasons. 
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